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1. Introduction 

In this contribution we report on the calculation of the next-to-leading 
order (NLO) DGLAP evolution kernels that was started in refs. [1,2]. The 
inclusive DGLAP splitting functions at NLO are known for a long time, they 
were first calculated in early 1980’s [3-6]. Today also the next-to-next-to- 
leading order (NNLO) corrections are known for over ten years [7,8] so why 
the interest in NLO splitting functions? The reason for this is the fact that 
all the above mentioned calculations were performed in the MS scheme and 
by their nature are inclusive. A direct motivation for the current work comes 
from Monte Carlo simulations of parton cascades in form of parton shower 
(PS) generators. In particular from the effort of the Cracow group that aims 
to include full NLO corrections in the parton shower simulations [9-12]. 
This means not only corrections to the hard matrix element, as is done in 
approaches like MC@NLO [13] and POWHEG [14], but also the full NLO 
corrections to the parton cascades. 

The big difficulty with including NLO corrections in parton shower al¬ 
gorithms is connected with the nature of parton shower generators. Specifi¬ 
cally in PS approach all the subsequent parton emissions need to be treated 
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in fully exclusive way, meaning that their 4-momenta are not integrated 
(including transverse degrees of freedom). For comparison in a standard in¬ 
clusive QCD calculations there are two elements, Wilson coefficient (parton 
level hard matrix element) and parton distribution functions (PDFs), that 
sum up all the above mentioned emissions in an inclusive way. Moreover if 
we go beyond LO in PS approach we need additional prescription to avoid 
double counting of contributions in hard matrix element and in the shower. 
This is referred to as the NLO matching and it was a major improvement for 
PS in the last decade. It was originally done with MC@NLO method and 
shortly after the POWHEG method was introduced. Recently an alterna¬ 
tive solution from the Cracow group is also available referred to as KrkMLO 
method [12]. This new solution is distinct from the earlier methods in a 
number of ways, one of them being the use of different factorization scheme. 
Instead of using traditional MS scheme that was designed for simplicity of 
inclusive calculations this method uses a dedicated factorization scheme, 
referred to as Monte Carlo (MC) scheme, which was designed especially for 
Monte Carlo simulations of parton cascades. The purpose of introducing 
this scheme was not only NLO matching but rather simplification of a more 
difficult task that is including NLO corrections in the shower and maybe 
even later matching of the NLO shower with NNLO matrix element. One 
of the necessary ingredients needed for this purpose are NLO evolution ker¬ 
nels calculated in the MC scheme. We first calculated the real emission 
corrections to the non-singlet NLO splitting functions [1] as these were re¬ 
quired for proof of concept of the methodology to include NLO corrections 
in the shower [9, 10]. Recently we complemented the scheme by providing 
prescriptions for computing virtual corrections and also calculated them for 
the non-singlet case [2]. Currently we are hnishing the calculation of the 
NLO splitting functions and in this contribution we are reporting on the 
progress. 


2. Monte Carlo friendly regularization scheme 

The regularization prescription for computing NLO splitting functions 
is one of the ingredients of the Monte Carlo factorization scheme [12]. It 
is based on the prescription of Curci, Furmnski and Petronzio (CFP) in¬ 
troduced in [5] during first NLO kernel calculation in x space. The most 
important features of the CFP scheme are (i) the use of axial gauge; and (ii) 
principal value (PV) regularization for the axial gluon propagators [leading 
to infra red (IR) singularities], namely: 

1 _ In 

^ In ^ [ln]pv ’ [ln]pv {In)'^ + {pi)' 
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(NLO-qq-d-real) (NLO-qq-d-virt) 


Fig. 1: Selected non-singlet diagrams contributing to the Pqq kernel. 


where n is the axial gauge vector, p is an external reference momentum and 5 
is the PV geometrical regulator. Unfortunately the use of PV prescription 
in this way leads to complicated patterns for cancellation of singularities 
between real and virtual diagrams, which is unacceptable if we want to use 
them for the construction of PS Monte Carlo. 

We illustrate it on the example of the non-singlet diagram NLO-qq-d 
displayed in Fig. 1. If calculated using CFP (MS ) prescription the virtual 
graph gives (we provide here only the singular terms); 




(2) 


and the corresponding real graph reads: 
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where Pqq{x) = is the LO qq kernel, and Iq and Ii represent single and 
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double logarithmic singularities regulated by the PV regulator 5: 


Io= [ dx^^ = -ln6 + 0{S), 

Jo mpv 

Jo [x\pv 2 24 


(4) 


We can see that higher order e poles cancel between the corresponding 
real and virtual graphs, which is probably impossible to implement in a MC 
shower program. 

Because of this fact we have modihed the CFP prescription by changing 
the way PV regularization is applied. Instead of applying it only to the gluon 
propagators we apply it to all the singularities in the plus variable, l~^ = 
nl/np, see [2] . This results in simplihcation of the real-virtual cancellations 
which can be seen by examining results for NLO-qq-d graphs in the new 
prescription (referred as NPV). The virtual graph is given by: 


-pvirt NPV 
^ NLO-qq-d 




Pqq{x) 

Pqq{x)( 


— 4/o — 2 ln(l — x) — ln(x) + ^ 

- 6/i -|- 6/o ln(l — x) -|- 2Iq ln(x) 
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and the real graph is; 


-preal NPV 
^ NLO-qq-d 


Pqq{x) 


^2/i + 4/o ~ 2/o lii(l 


x) + 2Ioln(x) 


( 6 ) 


One can easily check that upon adding real and virtual graphs in both 
prescriptions leads to the same final results. However, from the MC point 
of view the situation in the NPV prescription is hugely improved. There 
are no triple e poles, the double poles occur only in the virtual graph and 
the real one is free of them.^ 

Actually a general statement is valid. ^ In case of the NPV regularization: 

(i) poles in virtual and real graphs are absent (replaced by 1/e In^ d-type 
structures), and (nearly) all e^ poles originate from the virtual graphs alone, 

(ii) all real emission graphs (contributing to NLO splitting functions) feature 
only single e poles, (hi) sum of real and virtual graphs for the same topology 
reproduces the corresponding sum calculated using CFP prescription. 

^ More specifically, the 1/e^ terms have been replaced by l/el+ J-type terms, regulated 
simultaneously by t and 5. 

^ With some exceptions related to running of the coupling constant and CFP projection 
operators. 
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There are some additional subtleties connected to the graphs contribut¬ 
ing to the running of the strong coupling (e.g. NLO-gg-f), these are partly 
discussed in ref. [ 2 ] . 

Using NPV regularization prescription is the first step on the way to the 
results in the Monte Carlo scheme. For the MC shower simulations we need 
purely 4-dimensional calculations and as we can see above the presented 
results still feature single (and partly double) e poles. This can be however 
avoided by introducing additional cut-off regularization for the overall scale 
integration. It was already done in [1] when we calculated real emission 
contributions to the non-singlet splitting functions and won’t be discussed 
here. 


3. Status of the calculation 

Currently we are finishing the calculation of the singlet splitting func¬ 
tions using the prescription described in the previous section and introduced 
in [ 1 , 2 ]. At the moment we have recalculated the full Pgg and Pqg kernel on 
the inclusive level and we have obtained a perfect agreement with the liter¬ 
ature. Below we present the inclusive “parton density” (defined as in [15]) 
for the case of Pgg. It is obtained as a sum of contributions from individual 
diagrams in Fig. 2 : 


^99 = ^ (1^) IJ Pggi^) (y - 8^0 - 81n(l - x) - 41n(x)^ 
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where: 


Pgg{x) 

S2{x) 


(1 — X + 

x(l — x) ’ 



dz , 
— in 

2 : 
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The kernel is now easily extracted from Tgg by taking twice its residue (twice 
as we use n = 4 + 2e); additionally we exclude factor which is the 

usual convention, then: 
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This result reproduces the well known Pgg kernel, that can be found in 
ref. [16] or in the original Furmanski and Petronzio paper [6]. The literature 
provides the single e terms defining the splitting functions but some of the 
double pole terms from the inclusive densities can also be cross-checked with 
the results in [17]. 

Similarly for the case of Pqg the inclusive density is given by (the con- 
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tributing diagrams are listed in Fig. 3): 
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where Pqg{x) = x^ + (1 — x)^ is the corresponding LO kernel. The resulting 
splitting function is equal to: 
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which agrees with the literature [6,16]. 

The double e pole terms in the inclusive densities T originate from two 
sources. First is building up of the running coupling constant, the second is 
artificial and connected with the subtraction defined through the projection 
operator (P in CFP [5]), entering via countergraphs like NLO-qg-il or NLO- 
gg-i. In the inclusive calculations the higher e poles are not important as 
splitting functions do not depend on them (they are defined by the residue 
- single pole part), making this kind of projection operator perfectly well 
suited for the inclusive computation. However, if we want to use the unin¬ 
tegrated distributions (exclusive version of F) in the Monte Carlo program 
this will become a problem. We can solve it by using a dedicated Monte 
Carlo factorization scheme [9,12] which redefines projection operators, how¬ 
ever, this is beyond the scope of this work and we will not investigate it here. 
Let us just mention that similar idea has been investigated by Oliveira et 
al. [18,19] in the context of physical factorization scheme allowing for better 
convergence of perturbative series. 

Currently we are working to complete the calculation for Pgq and Pqq 
kernel. Diagrams contributing to these calculations can be seen in Figs. 4 
and 5. Computation for Pgq is analogical to the one for Pqg. In case of Pqq 
the first (Born) level contributions appear only at order 0(ag) making this 
calculation very easy. 


4. Summary 

We have reported on the progress of re-calculation of the NLO splitting 
function using new NPV regularization prescription [2]. We have already 
calculated the Pgg and Pqg (and earlier Pqq) NLO splitting functions and 
we have reproduced the know inclusive results from the literature. The 
calculations for the Pgq splitting function is already advanced and the partial 
results we have obtained so far are also in agreement with the literature. 

The inclusive results presented here are not our main interest and serve 
us mainly as a cross-check of our calculations. The main interest for us are 
the unintegrated distributions in the new scheme that will be better suited 
for the Monte Carlo simulations. These distributions will be presented else¬ 
where. 
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Fig. 2: List of diagrams contributing to the Pgg kernel. The graphs d,f,sl,dF 
and fF contribute in two versions, depending on the location of the cut. 


REFERENCES 


[1] S. Jadach, A. Kusina, M. Skrzypek, and M. Slawinska, “Two real parton 
contributions to non-singlet kernels for exclusive QCD DGLAP evolution,” 
JHEP 08 (2011) 012, 1102.5083. 





































10 


Epiphany2015 printed on July 6, 2015 



(NLO-qg-hl) - (NLO-qg-il) 





1 

r 1 U 


s ^ i 

(NLO-qg-d2) 

o ^ 

o ^ 

D |C 

(NLO-qg-f) 



(NLO-qg-e 


Fig. 3: List of diagrams contributing to the Pqg kernel. The graphs dl, d2 
and f contribute in two versions, depending on the location of the cut. 




(NLO-gq-b) 


(NLO-gq-e2) 


Fig. 4: List of diagrams contributing to the Pgq kernel. The graphs dl, d2 
and f contribute in two versions, depending on the location of the cut. 




[2] O. Gituliar, S. Jadach, A. Kusina, and M. Skrzypek, “On regularizing the 
infrared singularities in QCD NLO splitting functions with the new Principal 
































































Epiphany2015 printed on July 6, 2015 


11 



(NLO-qqb-hl) - (NLO-qqb-il) 



(NLO-qqb-b) 


Fig. 5: List of diagrams contributing to the Pqq kernel. 


Value prescription,” Phys.Lett. B732 (2014) 218-222, 1401.5087. 

[3] E. G. Floratos, D. A. Ross, and C. T. Sachrajda, “Higher order effects in 
asymptotically free gauge theories: the anomalous dimensions of Wilson 
operators,” Nud. Phys. B129 (1977) 66-88. 

[4] E. G. Eloratos, D. A. Ross, and C. T. Sachrajda, “Higher order effects in 
asymptotically free gauge theories. 2. flavor singlet Wilson operators and 
coefficient functions,” Nud. Phys. B152 (1979) 493. 

[5] G. Curd, W. Furmanski, and R. Petronzio, “Evolution of parton densities 
beyond leading order: the non-singlet case,” Nud. Phys. B175 (1980) 27. 

[6] W. Furmanski and R. Petronzio, “Singlet parton densities beyond leading 
order,” Phys. Lett. B97 (1980) 437. 

[7] S. Moch, J. A. M. Vermaseren, and A. Vogt, “The three-loop splitting 
functions in QCD: The non-singlet case,” Nud. Phys. B688 (2004) 101-134, 
hep-ph/0403192. 

[8] A. Vogt, S. Moch, and J. A. M. Vermaseren, “The three-loop splitting 
functions in QCD: The singlet case,” Nud. Phys. B691 (2004) 129-181, 
hep-ph/0404111. 

[9] S. Jadach, M. Skrzypek, A. Kusina, and M. Slawinska, “Exclusive Monte 
Carlo modelling of NLO DGLAP evolution,” PoS RADCOR2009 (2010) 
069, 1002.0010. 

[10] M. Skrzypek, S. Jadach, A. Kusina, W. Placzek, M. Slawinska, et al, “Fully 
NLO Parton Shower in QCD,” Acta Phys.Polon. B42 (2011) 2433-2443, 
1111.5368. 

[11] S. Jadach, A. Kusina, W. Placzek, M. Skrzypek, and M. Slawinska, 
“Inclusion of the QCD next-to-leading order corrections in the quark-gluon 
Monte Carlo shower,” Phys.Rev. D87 (2013) 034029, 1103.5015. 

[12] S. Jadach, W. Placzek, S. Sapeta, A. Siodmok, and M. Skrzypek, “Matching 
NLO QCD with parton shower in Monte Carlo scheme - the KrkNLO 
method,” 1503.06849. 

[13] S. Frixione and B. R. Webber, “Matching NLO QCD computations and 
parton shower simulations,” JHEP 06 (2002) 029, liep-ph/0204244. 












12 Epiphany2015 printed on July 6, 2015 


[14] P. Nason, “A new method for combining NLO QCD with shower Monte 
Carlo algorithms,” JHEP 11 (2004) 040, hep-ph/0409146. 

[15] O. Gituliar, “Higher-Order Corrections in QCD Evolution Equations and 
Tools for Their Calculation,” 1403.6897. 

[16] R. K. Ellis and W. Vogelsang, “The evolution of parton distributions beyond 
leading order: the singlet case,” hep-ph/9602356. 

[17] G. Heinrich, Improved techniques to calculate two-loop anomalous dimensions 
in QCD. PhD thesis, Swiss Federal Institute of Technology, Zurich, 1998. 

[18] E. Oliveira, A. Martin, and M. Ryskin, “Treatment of the infrared 
contribution: NLO QED evolution as a pedagogic example,” Eur.Phys.J. 
C73 (2013) 2534, 1305.6406. 

[19] E. Oliveira, A. Martin, and M. Ryskin, “Physical factorisation scheme for 
PDFs for non-inclusive applications,” JHEP 1311 (2013) 156, 1310.8289. 


